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LONG-TERM  GOALS 

To  quantify  and  understand  the  effects  of  aggregation  dynamics  on  the  distribution  of  particles  in  the 
bottom  boundary  layer,  and  to  understand  how  the  properties  of  particles  (composition,  shape,  and 
internal  structure)  affect  their  optical  and  acoustical  properties. 

OBJECTIVES 

•  Obtain  direct  measurements  of  the  turbulent  Reynolds  stress  within  the  wave  boundary  layer 
(WBL) 

•  Obtain  concurrent  velocity,  turbulence  and  seafloor  microtopography  (bedforms) 
measurements  to  constrain  the  fluid  dynamical  environment  within  which  the  particle  size 
distribution  evolves. 

•  Improve  state-of-the-technology  for  measuring  boundary  layer  velocity  and  turbulence  profiles 
(Pulse  Coherent  Doppler  Profiler)  and  seafloor  bedform  measurements  (Rotary  Sidescan  and 
Pencil-beam  sonars) 

APPROACH 

Our  project  is  focused  on  obtaining  direct  measurements  of  the  turbulent  Reynolds  stress  within  the 
wave  boundary  layer  (WBL),  which  has  not  been  possible  in  previous  field  studies.  The  measurements 
of  stress  are  combined  with  detailed  topographic  maps  of  the  bedforms  to  understand  the  relations 
between  the  bedforms  and  the  turbulent  flow  over  the  bedforms.  This  is  motivated  by  the  order-of- 
magnitude  difference  between  relatively  small  stresses  that  have  been  measured  above  the  WBL  and 
the  much  larger  stresses  that  have  been  inferred  within  the  WBL.  The  discrepancy  between  measured 
and  modeled  stresses  above  the  WBL  has  been  attributed  to  processes  within  the  WBL.  Hill  et  ah, 
(2001)  inferred  the  large  stresses  within  the  WBL  play  a  dominant  role  in  controlling  particle  size  and 
thus  strongly  influence  particle  settling  velocities,  concentrations,  and  optical  properties. 

The  measurements  capitalized  on  multi-frequency  pulse-coherent  Doppler  profiling  (PCDP)  sonars 
developed  recently  at  WHOI  by  Peter  Traykovski  and  Fred  Jaffre  (2011)  based  on  a  high  frequency 
sonar  board  designed  by  Gene  Terray  and  Tom  Austin  and  Peter  Traykovski.  The  multi-frequency 
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methodology  (Hay  and  Zedel,  2008)  removes  the  ambiguity  associated  with  conventional  pulse- 
coherent  sonar  measurements.  The  combination  of  range-gated  vertical  (transmit-receive)  and  slanted 
(receive-only)  acoustical  beams  produces  vertically  resolved  measurements  of  the  horizontal  and 
vertical  velocity  throughout  the  wave  boundary  layer  and  the  lower  part  of  the  overlying  current 
boundary  layer  (Figure  1).  The  effective  sample  volume  is  the  product  of  the  beam  patterns  along  the 
acoustical  transmission  and  reception  paths,  and  is  small  because  of  the  narrow  vertical  beam.  The 
new  sensors  produce  high-precision  velocity  measurements  along  with  order-one-meter  range,  sub¬ 
centimeter  spatial  resolution,  and  a  sample  rate  of  approximately  18  Hz.  With  proper  configuration, 
these  sensor  characteristics  are  sufficient  for  measurements  of  both  the  oscillatory  velocity  and  the 
turbulent  velocity  fluctuations  in  the  WBL  at  field  scales. 

WORK  COMPLETED 

The  instrumentation  was  deployed  in  1 1  m  water  depth  at  the  Martha’s  Vineyard  Coastal  Observatory 
(MVCO)  in  September  and  October  of  201 1  for  a  period  of  one  month  (Figure  2).  Data  processing  and 
analysis  was  conducted  and  the  results  were  presented  at  the  2012  American  Geophysical  Union  Ocean 
Sciences  Meeting  in  poster  entitled  “Observations  of  Wave  Boundary  Layer  Turbulence  and  Sediment 
Suspension”  The  2011  ONR  annual  report  describes  the  details  of  the  equipment  development  and 
preparation  for  this  work  (http://www.onr.navy.mil/reports/FYl  l/lotrayk2.pdf). 

RESULTS 

The  rotary  sonar  imagery  shows  orbital  scale  bedforms  in  the  coarse  sand  (D50=  65  0  pm)  present  in  the 
rippled  scour  depressions  at  MVCO.  These  bedforms  have  wavelengths  that  vary  from  0.40  to  1.2  m 
and  scales  with  the  orbital  diameter  of  the  waves.  Figure  3  shows  an  example  with  0.7  m  wavelength. 

Horizontal  and  vertical  velocity  profiles  from  the  PCDP  show  energetic  wave  flows  and  bursts  of 
turbulence  shed  from  the  ripple  crests  (Figure  4).  Inertial  dissipation  method  fits  to  spectra  from  both 
time  series  analysis  and  direct  spatial  analysis  reveal  a  relatively  low  friction  factor  (fe  =  0.01)  and 
consistent  with  grain  roughness  in  the  well-known  Swart’s  (1976)  friction  factor  (Figure  5). 

A  two  dimensional  analysis  of  the  turbulent  flow  over  the  ripples  estimated  by  using  the  ripple 
migration  past  the  PCDP  as  a  proxy  for  position  along  the  ripple  crest  allows  analysis  of  the  strength  of 
the  eddies  shed  from  the  ripple  crests  (Figure  6).  These  eddies  produce  stresses  that  are  highly  out  of 
phase  with  free  stream  velocity  and  are  thus  high  dissipative.  An  analysis  of  the  dissipation  due  to  the 
total  ripple  roughness,  including  the  eddies  produces  a  friction  factor  of  (fe  =  0.1)  consistent  with 
laboratory  results  over  orbital  scale  ripples.  This  represents  one  of  the  first  field  measurements  of 
wave  boundary  layer  turbulent  flow  over  orbital  scale  ripples. 

IMPACT/APPLICATIONS 

Operational  seagoing  systems  often  depend  on  optical  and  acoustical  properties  of  suspended  particles 
in  the  water  column.  Understanding  the  processes  that  regulate  the  particle  characteristics  and 
understanding  the  optical  and  acoustical  signatures  of  suspended  particles  are  essential  in  order  to 
predict  the  performance  of  these  operational  systems.  The  development  of  instrumentation  that  can 
resolve  wave  boundary  layer  turbulence  has  potential  to  make  significant  improvements  to  our 
understanding  of  boundary  layer  and  sediment  processes. 
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RELATED  PROJECTS 


Much  of  the  development  work  and  equipment  purchases  described  in  this  work  was  funded  by  a 
DURIP  proposal  entitled  “Instrumentation  for  Measuring  Nearshore  Morphologic  Change  and 
Hydrodynamic  Forcing”  (00014-10-10768)  and  instrumentation  will  also  be  used  in  the  “Dynamics  of 
sandwaves  under  combined  wave-current  forcing  and  mine  burial  processes”  (N000 14-1 1-10291)  and 
“Multi-Scale  (cm  to  km)  Hydrodynamic  and  Morphologic  Interactions  in  Tidal  Inlets”  (N000 14-10- 
10376) 


Pulse  Coherent 
Doppler  Profiler 


Figure  1.  Schematic  diagram  of  the  pulse-coherent  Doppler  profiler  developed  at  WHOI .  The  red 
shaded  area  represents  the  region  of  1  cm  resolution  profiles  of  u,  v  and  w  at  turbulence  resolving 
rates  (18  Hz).  An  additional  pair  of  angled  transducers  is  located  in  the  plane 

perpendicular  to  the  pair  shown. 
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Figure  2.  Photograph  of  instrumentation  just  before  deployment  at  MVCO  in  September  2011, 
showing  Pulse  Coherent  Doppler  Profiler  (PCDP),  Imagenex  rotary  sidescan  sonar,  Imagenex  2- 
axis  rotary  pencil  beam  sonar,  Aquatec  ABS  and  Nortek  Vector  AD  Vs. 
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Rotary  Sonar  Date:17-Sep-201 1  23:17:20 
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Figure  3.  Rotary  Sidescan  Image  showing  70  cm  wavelength  orbital  scale  (l  ~  0. 75  ( 2A ))  ripples, 

where  2A  is  the  orbital  diameter 
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Figure  4.  left,  top  to  bottom:  vertical  velocity  profiles,  horizontal  velocity  profiles  and  pressure 
right:  vertical  velocity  spectra  at  (blue)  0.45  in,  (red)  0.05  m  above  the  bed  and  (magenta)  spectra 

based  on  direct  wavenumber  spatial  spectra. 
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Vertical  Profile  of  s 


s  calculated  from  Inertial  Dissipation  Method  (IDM) 
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Time  series  of  D  =  p  f  e  dz 


Figure  5.  left:  time  averaged  profdes  of  dissipation  estimated  using  the  inertial  dissipation  method. 

right:  time  series  of  dissipation  profiles.  Bottom:  time  series  of  depth  integrated  dissipation 
Dissipation  associated  with  measured  turbulent  fluctuations  is  low  (f.  =  0.01)  and  consistent  with 

Grain  Roughness  in  Swart’s  Friction  Factor 
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Vorticity  and  Velocity  Deficit  using  ripple  migration  as  proxy  for  x 
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Vorticity  and  Velocity  Deficit  using  ripple  migration  as  proxy  for  x 


Figure  6.  Left  panel:  phase  averaged  flow  field  and  vorticity  at  a  phase  of 330°  past  (or  30°  before) 
maximum  onshore  velocity.  Positive  vorticity  is  visible  in  the  lee  of  ripple.  This  vorticity  was 
advected  from  the  separation  in  the  period  between  180  and  270°.  Right  panel:  phase  averaged  flow 
field  and  vorticity  at  a  phase  of  60° past  maximum  onshore  velocity  showing  negative  vorticity 
consistent  with  the  separation  of  flow  in  the  lee  of  the  ripple. 
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Figure  7.  Dissipation  associated  with  measured  turbulent  fluctuations  is  low 
(fe  =  0.01)  and  consistent  with  grain  roughness  in  Swart’s  Friction  Factor  (magenta  dots).  When 
this  dissipation  associated  with  measured  turbulent  fluctuations  is  placed  on  the  friction  factor 
graph  at  a  roughness  associated  with  ripple  height,  the  data  falls  well  below  Swart’s  formulation 
(magenta  circles).  Using  a  momentum  deficit  approach  in  the  near  bed  region  produces  a  slightly 
higher  friction  factor  (red  triangle,  fe  =  0.03),  but  in  order  to  produce  a  friction  factor  consistent 
with  Swart  s  formulation  the  dissipation  from  the  large  eddies  shed  from  the  ripple  crests  must  be 

included  (blue  triangle,  fe  =  0.1). 
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